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Newly born massive magnetars are generally considered to be produced by binary neutron star (NS) mergers, 
which could give rise to short gamma-ray bursts (SGRBs). The strong magnetic fields and fast rotation of 
these magnetars make them promising sources for gravitational wave (GW) detection using ground based GW 
interferometers. Based on the observed masses of Galactic NS-NS binaries, by assuming different equations of 
state (EOSs) of dense matter, we investigate the stochastic gravitational wave background (SGWB) produced 
by an ensemble of newly born massive magnetars. The massive magnetar formation rate is estimated through: 
(i) the SGRB formation rate (hereafter entitled as MFR1); (ii) the NS-NS merger rate (hereafter entitled as 
MFR2). We find that for massive magnetars with masses Mmg = 2.4743Mo, if EOS CDDM2 is assumed, the 
resultant SGWBs may be detected by the future Einstein Telescope (ET) even for MFR1 with minimal local 
formation rate, and for MFR2 with a local merger rate 62(0) < 10 Mpc?Myr !. However, if EOS BSk21 is 
assumed, the SGWB may be detectable by the ET for MFR1 with the maximal local formation rate. Moreover, 
the background spectra show cutoffs at about 350 Hz in the case of EOS BSk21, and at 124 Hz for CDDM2, 
respectively. We suggest that if the cutoff at ~ 100 Hz in the background spectrum from massive magnetars 
could be detected, then the quark star EOS CDDM2 seems to be favorable. Moreover, the EOSs, which present 


relatively small TOV maximum masses, would be excluded. 


PACS numbers: 04.30.-w, 97.60.Jd, 26.60.Kp, 04.30.Db 


I. INTRODUCTION 


Short gamma-ray bursts (SGRBs) are generally considered 
to be arising from the coalescence of either a neutron star- 
neutron star (NS-NS) binary or a neutron star-black hole (NS- 
BH) binary (see for a recent review). Mergers of these 
compact binaries produce strong gravitational wave (GW) 
emissions, making them promising sources for GW detection 
using the ground based GW interferometers such as LIGO, 
VIRGO, GEO600, advanced LIGO (aLIGO), and the future 
Einstein Telescope (ET) (2, B). The merger product of a NS- 
BH binary is of course a stellar-mass BH. On the other hand, 
the remnant of a NS-NS merger is still an open question. De- 
pending on the total mass of the binary system and the NS 
equation of state (EOS), the NS-NS merger product may be 
either of the following four possibilities : (1) a stellar- 
mass BH; (2) a differential rotation supported unstable hyper- 
massive NS, which will collapse into a BH in a few tens mil- 
liseconds; (3) a centrifugal force supported temporarily stable 
massive NS, which will collapse into a BH when the NS is 
spun down; (4) an eternally stable massive NS. The massive 
NS remnants are suggested to possess strong surface dipole 
and internal toroidal magnetic fields, which are amplified due 
to various mechanisms, such as Kelvin-Helmholtz instability 

, magnetorotational instability [9] during/after the merger, 
qa — € dynamo in the nascent millisecond NS [10], and the 
combined effect of r-mode and Tayler instabilities [11]. Ob- 
servationally, the existence of extended emissions (12) , X-ray 


flares [13]. and internal x-ray plateaus in a large sample 
of SGRBs x-ray lightcurves support the idea that the central 
object of some SGRBs could be a highly magnetized, mil- 
lisecond rotating NS. 


Strong GW emission is expected in the final inspiral pro- 
cess of a NS-NS binary. Moreover, if the merger product is 
either (3) or (4), the remnant can also produce strong long- 
lasting GW signals, though their strengths may be relatively 
weak. Generally, the fast rotating, massive NS can emit 
GWs because of nonaxisymmetric instabilities, e.g., dynam- 
ical bar-mode instability [15], r-mode instability , f-mode 
instability (17). On the other hand, strong internal magnetic 
fields of the massive magnetar can lead to nonaxisymmet- 
ric quadrupole deformation, which could also produce GW 
emission fig orn. The amplitude of the magnetically induced 
GW signal is proportional to the quadrupole ellipticity [18]. 
which mainly depends on the EOS, the magnetic energy, and 
the interior magnetic field configuration of the NS (see, e.g. 
[18.20 2129). 

Superposition of the magnetically induced GW emissions 
from an ensemble of magnetars throughout the Universe can 
contribute to the astrophysical stochastic GW background 
(SGWB). The SGWB from the magnetic deformation of 
newly born magnetars has been discussed in many literature 
references (5d 3. However, in these papers, the magnetar 
mass is assumed to be a canonical value of 1.4Mo, which 
means these magnetars are eternally stable. Actually, the mag- 
netars produced by NS-NS mergers apparently have masses 


much larger than 1.4Mo, and they may not be always stable 
[e.g., product (3) mentioned above]. For a single source, the 
collapse of the magnetar will lead to the cease of GW emis- 
sion at a certain frequency. Hence, in order to derive a real- 
istic SGWB produced by the magnetic deformation of newly 
born massive magnetars, both the EOS of dense matter and 
the masses of magnetars should be taken into account. In this 
paper, we reconsider the SGWB produced by magnetic defor- 
mation of the newly born, massive magnetars based on typi- 
cal NS and quark star (QS) EOSs, and the observed masses 
of the Galactic NS-NS binaries. The motivation of consider- 
ing the QS EOS is the recent statistical analysis of internal 
x-ray plateaus of SGRBs that shows that QS remnants might 
be more preferred than NSs BA. The paper is organized as 
follows: in Sec. II, we show how the GW signal from a single 
newly born massive magnetar is affected by the EOS. In Sec. 
III, we estimate the massive magnetar formation rate (MMFR) 
based on two different results: (1) the SGRB rate at redshift z 
suggested in B5]; (ii) the NS-NS merger rate as predicted by 
Regimbau and Hughes Bal. Results for the SGWBs produced 
by an ensemble of newly born massive magnetars are shown 
in Sect. IV. Conclusion and discussions are presented in Sect. 
V. 


Il. GW EMISSION FROM THE NEWLY BORN MASSIVE 
MAGNETAR 


We assume that all SGRBs are produced by NS-NS merg- 
ers, and the merger remnants are either temporarily or eter- 
nally stable massive NSs/QSs. Generally, in the merger pro- 
cess, only € 10-?Me materials are ejected from the NS-NS 
binary system (371. Hence, the total rest mass of the sys- 
tem is basically conserved, e.g., Mim = M; + Mi», where 
M; represents the rest mass of the massive magnetar rem- 
nant, M, ; and M; are the rest masses of the two NSs, respec- 
tively. Based on the approximate relation between rest 
and gravitational masses, by assuming that the extragalactic 
NS-NS binaries have the same gravitational mass distribution 
as the Galactic NS-NS binary population, one can easily esti- 
mate the distribution of gravitational mass Mmg for the mas- 
sive magnetar remnants i Bá. 

Until now there are six Galactic NS-NS binaries that have 
a relatively accurate measured gravitational mass for each NS 
in the binary system; they are PSR J0737-3039 (with gravita- 
tional masses My = 1.3381Me and M,» = 1.2489 Me for the 
two NSs, respectively), PSR B1534+12 (Mg; = 1.3332Mo 
and M,» = 1.3452Mo,), PSR B1913+16 (M,; = 1.4398Mo 
and M,» = 1.3886M,), PSR B2127+11C (Mg; = 1.358Mo 
and Ma2 = 1.354M 5), PSR J1906+0746 (Mg, = 1.248Mo 
and Mg = 1.365Mo), and PSR J1756-2251 (Mg; = 1.40Mo 
and M,» = 1.18Mo) . Therefore, the gravitational mass of 
the massive magnetar remnant is 2.4046Mo as inferred from 
the binary system PSR J0737-3039, 2.4845 Mo inferred from 
PSR B15344+12, 2.6154Mg inferred from PSR B1913+16, 
2.5139Mo, inferred from PSR B2127+11C, 2.4276Mo in- 
ferred from PSR J1906+0746, 2.3996M, inferred from PSR 
J1756-2251. The average gravitational mass of the massive 


magnetars is thus Mmg = 2.4743Mo, and we take this value as 
the typical mass for the massive magnetars in the remainder 
of this paper. 

For a nonrotating NS/QS, the maximum gravitational mass 
that it could sustain is the Tolman-Oppenheimer-Volkoff 
(TOV) maximum mass, Mroy, which is determined by the 
NS/QS EOSs. However, centrifugal forces due to the uniform 
rotation of the merger remnant could increase the maximum 
sustainable gravitational mass. Li et al. 13.4] calculated equi- 
librium sequences of uniformly rotating NS/QS configura- 
tions with a spin frequency increasing from 0 to the Keplerian 
spin limit and obtained analytical expressions for the max- 
imum gravitational mass Mg max, the corresponding equilib- 
rium radius Reg (in kilometers), and the corresponding max- 
imum moment of inertia Imax of a NS/QS with a spin period 
P (in milliseconds), which, respectively, have the following 
form: 


Mg,max = Mrov(1 + a PP); (1) 
Req = Ct APP; (2) 
Imax = M, R? : (3) 


g,max‘Yeq T 4 e kP)’ 


where M, max and Myoy are measured in solar masses. The 
fitting parameters a, B, A, B, C, a, q, and k are EOS depen- 
dent. For the typical NS (BSk21 [40]) and QS (CDDM2 [41]]) 
EOSs considered in this paper, the specific values of these pa- 

rameters as well as Mroy can be found in Table 1 of [34]. 
From Eq. (1), one can define the collapse frequency, veon = 
1/Pcou, below which the massive magnetar with a gravita- 
tional mass Mmg = Mg max(Pcon) will immediately collapse 
into a BH. The specific form of veon is [5H7] 
| æœMrovy | (4) 

Ycoll Mmg = Mroy . 


If vecon € 0 G.e., Mmg < Mroy), the massive magnetar is eter- 
nally stable. However, if 0 < Yeon < vi (with vj represents 
the initial spin frequency), the massive magnetar is temporar- 
ily stable, and it will collapse into a BH when the star spins 
down to Veo due to GW emission and magnetic dipole radi- 
ation (MDR). Lastly, if veon > vi, the massive magnetar will 
collapse into a BH immediately after it is born. Remarkable 
GW emissions from the central remnant are expected only 
in the first two cases (i.e., eternally stable and temporarily 
stable magnetars). For the EOSs BSk21 and CDDM2 con- 
sidered here, the massive magnetar remnant with an initial 
spin at the Keplerian limit should be temporarily stable be- 
cause 0 < Yeo < vi (see below). Other EOSs that provide 
Mqiov > 2.4743Mo will result in an eternally stable massive 
magnetar, and continuous GW emission extended to lower fre- 
quencies. 

The newly born massive magnetar spins down mainly 
through MDR and magnetically induced GW emission. 
Therefore, the evolution formula for the angular frequency w 
of the magnetar can be written as 


BRs% — 32Geglo? 
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(5) 


w= 


where Bg is the surface dipole magnetic field at the mag- 
netic pole, R the radius, and J the moment of inertia of 
star. By adopting different interior magnetic field configura- 
tions and stellar interior structures, the magnetically induced 
quadrupole ellipticity eg has been calculated in many litera- 
ture references (e.g., fig ). Some nonlinear numer- 
ical simulations show that the interior magnetic field prob- 
ably has a poloidal-toroidal twisted-torus shape (42). How- 
ever, even for this configuration, the dominated one is usually 
the toroidal field component. In the toroidal-dominated case, 
€g is related to the volume-averaged strength of the toroidal 
field B, (23), which is hard to be determined directly. Gen- 
erally, B,/Ba ~ 5-100 with Bg the dipole magnetic field of 
the magnetar is proposed following the observations of gi- 
ant flare from SGR 1806-20 [19], free precession of magne- 
tar 4U 0142+61 (43h, x-ray afterglows of some SGRBs 
and lightcurves of superluminous supernovae 43]. For a NS 
the ellipticity can be estimated as eg ~ 1074(B,/10!® Gy 
46]. While for a QS, if it is in the two-flavor color su- 
perconductivity a Dl the ellipticity is approximated as 
€3 ~ 7 x 10-*(B,/10!° G) for the mass and radius adopted 
thereinafter [50]. On the other hand, eg can be constrained 
via analyzing the internal x-ray plateau afterglows of SGRBs 
[34]. Specifically, depending on the EOSs, the elliptic- 
ity of the massive magnetar is confined to be ey = 0.002 
for NS EOS BSk21 and eg = 0.004—0.007 for QS EOS 
CDDM2 [34]. Thereinafter, the representative ellipticities 
€g = 0.005 (the value with the best Kolmogorov-Smirnov 
test), and eg = 0.002 will be taken for EOSs CDDM2, and 
BSk21, respectively, while calculating the GW signal emitted 
by a single magnetar and the SGWB from the massive magne- 
tar population . The corresponding strength of the toroidal 
field is thus B, ~ 4.5 x 10!° (7.1 x 10/6) G for a NS (QS). 

The GW energy spectrum emitted by a single newly born 
massive magnetar can be estimated as 


E) 


dEgw | 32nG - 
Te earl, 


where ve = w/7 is the GW frequency at the source frame. One 
can also obtain the characteristic amplitude of the emitted GW 


as follows [51152]: 


ny ye 2 (7) 
C Ye = $$ 

vVdy./dt 
Ax  Gleg Y? 


where h(t) = ——77— is the GW strain amplitude, d is the dis- 
tance to the source. To assess the detectability of the GW sig- 
nal, we calculate the optimal (matched-filter) signal-to-noise 
ratio (SNR) as 


Ye. max h2 1/2 
S/N = | f Ao | (8) 
Vena vS n(Ve) 


! The rotation and temperature observations of pulsars disfavor the color- 


flavor-locked QS model (a8 [49]. 


where Ve min(=2Vcou) and Vemax(=2v;) are, respectively, the 
minimum and maximum GW frequencies emitted by the mag- 
netar, S (ve) is the one-sided noise power spectral density of 
the detector. The analytical expressions of S (ve) can be found 
in [53]] for aLIGO and ET. 


We do not follow instantaneous variations of the gravita- 
tional mass, radius, and moment of inertia with the spin-down 
of the massive magnetar, though all these quantities should 
actually decrease?. For simplicity, we take a typical gravita- 
tional mass Mmg, radius R, and moment of inertia J for the 
massive magnetar and assume they do not evolve with time 
during spin-down. For a magnetar with Mmg = 2.4743M,, its 
R and I are EOS dependent, which can be estimated as fol- 
lows. The radius R is approximately estimated by substituting 
the derived collapse period Poon = 1/Veon into Eq. (2). Then 
with R and Pon, using Eq. (3), the moment of inertia J can be 
obtained approximately. The resultant radius and moment of 
inertia of the 2.4743 Mg magnetar are, respectively, R = 12.66 
(16.31) km and J = 3.68 x 10? (5.50 x 109) g cm? if EOS 
BSk21 (CDDM2) is assumed. Obviously, R and J are under- 
estimated for the 2.4743 Mg magnetar that initially spins at the 
Keplerian limit Px. As a rough estimation, assuming a con- 
stant mass Mmg = 2.4743Mo, the ratio between the magnetar 
radii obtained at Px and at Poon is R(PxK)/R(Peon) « 1.1 (1.4) 
for EOS BSk20 (CDDM I1) (see Fig. 1 of ). Furthermore, 
with the spin-down of the magnetar, R(P)/R(Pcon) should de- 
crease and become equal to 1 when Peon is reached, where 
R(P) denotes the instantaneous radius of the magnetar with 
a spin period P. For EOSs BSk21 and CDDM2 considered, 
R(Pk)/R(Pcon) are not expected to vary too much from the 
above values. Following Eqs. (7) and (6), we have A. œ R, and 
dEgw dv. œ R? during the early period of spin-down when 
the GW emission is dominant. Consequently, our choice of a 
constant R will at most underestimate the characteristic ampli- 
tude h., and the background emission Qgw [see Eq. (14)] by 
a factor of 1.4 and 2, respectively. Hence, it is reasonable to 
take a constant R and / for a specific EOS in the calculations 
below?. 


Following Li et al. [34]. the dipole magnetic field of the 
massive magnetar is taken to be By = 10! G, and the ini- 
tial angular frequency is taken as wi = 27/Px. The values of 
Px for EOSs BSk21 and CDDM2 can be found in [34]. As- 
suming EOSs BSk21 and CDDM2, we show the characteristic 
amplitude A. of the GW signal versus the emitted frequency 
Ve in Fig.|1| The distance to the source is taken to be d = 100 
Mpc. The massive magnetar with Mmg = 2.4743Mg is tem- 
porarily stable for EOSs BSk21 and CDDM2, and it will not 
collapse until vco is reached. The collapse of the magnetar is 
manifested as a catastrophic cutoff in the emitted GW signal 
at 2v;gj, which is about 2453 Hz for EOS BSk21, and 868 


? During the spin-down process of a constant baryon mass massive magnetar, 
its gravitational mass decreases more slightly, in contrast to the radius and 
moment of inertia, which show very obvious decreases (see Fig. 1 of BA. 

3 The changes in Mmg, R, and / during the spin-down of massive NSs are also 
neglected in idm, since these effects are unlikely to significantly affect the 
evolutions of massive NSs and further their final results. 
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FIG. 1: The GW characteristic amplitude A, versus the emitted fre- 
quency ve, calculated by assuming NS EOS BSk21 (red line) and QS 
EOS CDDM2 (green line). For comparison, the rms strain noises for 
LIGO (solid line), VIRGO (dashed line), aLIGO (dotted line), and 
the future ET (dash-dotted line) are also presented [53]. 


Hz for EOS CDDM2, as shown in Fig. For EOS BSk21, 
the GW signal emitted by the massive magnetar extends from 
about 3322 Hz down to 2453 Hz. While for EOS CDDM2, 
the emitted GW is at lower frequency band, which covers the 
range from 1778 to 868 Hz. 


Since the strength of B, of a newly born magnetar is highly 
uncertain, in order to comprehensively show how B, could af- 
fect the detectability of GW signal from a single source, in 
Fig. the SNR S/N is plotted as a function of B,, whose 
range is ~ 5-100B, as discussed above. Obviously, with 
the increase of B, the SNR is gradually enhanced. However, 
for different EOSs, the evolution behaviors of S/N with B; 
differ significantly. Compared with EOS CDDM2, as B, in- 
creases, S/N shows a more obvious trend of getting saturated 
when EOS BSk21 is assumed. This is because for NS EOS 
BSk21, the ellipticity is more sensitive to the increase of B, 
(eg œ B? versus eg œ B, for QS EOS CDDM2). When eg 
is large enough, GW emission will dominate the spin-down; 
thus, we have S/N œ h(t)/ Ndv./dt = const. Moreover, us- 
ing the same detector, the derived SNR is higher for EOS 
CDDM2, since the sensitivities of the detectors are better at 
a relatively low frequency band. Assuming EOS BSk21, the 
SNRs of the GW emitted by the magnetar with a representa- 
tive ellipticity eg — 0.002 are 4.18 for ET (red filled star in 
Fig.|2) and 0.17 for aLIGO (red hollow star). While assuming 
EOS CDDM2 and eg = 0.005, we have S/N=16.65 for ET 
(green filled star in Fig. and S/N=0.72 for aLIGO (green 
hollow star). Adopting a single-detector search, the detection 
threshold is roughly S/N=8 . Hence, using ET, the emitted 
GW by the 2.4743 Mo. magnetar at 100 Mpc is undetectable if 


^ In newly born magnetars, B, ~ 10!7 G is possible (see, e.g. Eleda, 
since this strength is still lower than the virial limit by about an order of 


magnitude (26). 
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FIG. 2: The SNR S/N of the GW signal emitted by a single massive 
magnetar versus the magnetar's interior toroidal magnetic field B,. 
The SNRs are calculated by assuming NS EOS BSK21 (red lines), 
and QS EOS CDDM2 (green lines), respectively. The solid lines 
represent the SNRs obtained using future ET, while the dashed lines 
show the results derived using aLIGO. The dotted line indicates the 
detection threshold for a single-detector search. The stars show the 
SNRs calculated by adopting various representative ellipticities and 
for different detectors (see the text). 


EOS BSk21 is assumed even for B, = 10!’ G. For compari- 
son, when EOS CDDM2 is assumed, a detectable GW signal 
is expected if B, = 2.3 x 10!6 G (see Fig. [2), which may easily 
be achieved for newly born magnetars. Consequently, if future 
ET could detect the GW emitted by the 2.4743M. magnetar, 
then QS EOS CDDM2 will be more preferred. Moreover, ob- 
servation of the cutoff in the GW signal using ET may provide 
us an important channel to distinguish different EOSs. 


IIl. THE MASSIVE MAGNETAR FORMATION RATE 


As à quite rough estimation, the MMFR can be considered 
to be equal to the formation rate of SGRBs because we have 
assumed that only NS-NS mergers produce SGRBs and the 
merger products can only be temporarily stable or eternally 
stable massive magnetars. Obviously, this assumption will 
lead to an overestimation of the MMFR. Using the spectral 
peak energy-peak luminosity correlation for SGRBs, Yone- 
toku et al. [35] determined the redshifts of 72 BATSE SGRBs 
and obtained the relation between the SGRB formation rate 
and the redshift, which has the following form: 


PscgE(0X1-- 25$ , 0«z«0.67 


x ©) 
pscrp(O) x 1.676 , z 2 0.67 


PscRB(Z) = 


where Pscrp(0) is the local SGRB formation rate?. Hereafter, 
we refer to this as magnetar formation rate 1 (MFR1) and as- 
sume that the MMFR can be described by Eq. (9) up to z. ~ 6. 
The minimum SGRB formation rate at z = 0 is estimated to be 
PSGRB,min(0) = 1.15 x 107? events Mpc?yr^! by involving the 
geometrical correction of beaming angles B5). On the other 
hand, using the peak fluxes of 14 Swift SGRBs, redshifts, and 
beaming angles inferred from x-ray observations, Coward et 
al. obtained an upper limit for the local formation rate 
aS PSGRB,max(0) = 1.1 x 1076 events Mpc^?yr^! in the case of 
beamed emission. 

To estimate the MMFR, one can also equivalently esti- 
mate the NS-NS merger rate under the assumption that NS- 
NS mergers can only produce temporarily stable or eternally 
stable massive magnetars. Hereafter, we refer to the MMFR 
derived in this way as magnetar formation rate 2 (MFR2). As- 
suming that the NS-NS merger rate tracks the cosmic star for- 
mation rate (CSFR) with the time delay tg from formation of 
the NS binary to the final merger, the observed NS-NS merger 
rate at redshift z can be written as 


Pc) 
Ps (0) 


where ?(0) is the observed local merger rate per unit volume, 
whose value can be extrapolated by multiplying the Galactic 
NS-NS merger rate with the density of Milky Way-like galax- 
ies. Following Regimbau and Hughes Ba, we take two rep- 
resentative values for the local merger rate: (i) p?(0) = 10 
Mpc^?Myr^!, which is the upper limit of the local merger rate; 
(11) 68(0) = 1 Mpc ?Myr !, which represents the most proba- 
ble value for the rate. The NS-NS merger rate is related to the 
CSFR by the quantity b+ (z), which can be given as 


Pe(Z) = pe (0) x 


(10) 


pO. | E Pandu an 
+ Zf 

where o. is the CSFR. z and zp are the redshifts at which 
the NS-NS binary mergers and its progenitor binary initially 
formed, respectively. tą is the time difference between the 
formation of the progenitor binary and the compact binary, 
plus the merging time of the binary. It also represents the 
lookback time between z and zr, which has an approximate 
form fa = sa / Onl +z + QA— Onl zr)? + Qa). 
In this paper, the ACDM cosmological model is taken with 
the Hubble constant Hp = 70 km s !Mpc |, Qi = 0.3, and 
Q4 = 0.7. Based on the result of population synthesis (see 
(ad! and references therein), the probability distribution P(fa) 
for the time delay fg is given by 


P(ta) œ 1/ta with ta > To, (12) 


for some minimal delay time ro. For NS-NS binary, the min- 
imal delay time is assumed to be Tọ = 20 Myr, which corre- 
sponds to the evolution time from massive binaries to NS-NS 
binaries Bal. 


5 Since we use this formula as a rough estimation of the MMFR, the error 
bars in the exponent of (1 + z) and the expressions for pscna(0) are all 
neglected. 


For the CSFR, we use the result suggested in Hopkins and 
Beacom Bal. Based on the new measurements of the galaxy 
luminosity function in the UV and far-infrared wavelengths, 
they refined the previous models up to z. ~ 6 and obtained a 
parametric fit formula for the CSFR, which takes the follow- 
ing form [5g]: 


0.017 + 0.13z 


-1 -3 
Teas (3.333 Mayr Mpc ~~, (13) 


Èz) =h 


where h = 0.7. 


IV. RESULTS 


The SGWB is generally represented by the dimensionless 
quantity, Ocw(vos,), which describes the distribution of the 
GW energy density versus the GW frequency in the observer 
frame vobs. The SGWB produced by the magnetic deforma- 
tion of an ensemble of newly born massive magnetars is given 


by 


8zGvobs “upp 
Q obs) = : 
Gw( b: ) 3Hic f 


pMrR(Z) | dÉGw z 
(1+ ZE(Q,z) dve ” 


(14) 


Zlow 


where vos, = Ve/(1 + z), E(Q,z) = VOQn + zP + Qa, and 
Pmer(Z) is the MMFR, which can be substituted by Eqs. (9) 
or (10). Zupp and the Ziow are the upper and lower limits of the 
redshift integration, respectively. Zupp depends on the maximal 
redshift z. of the MMFR model and the maximal value of ve, 
1.€., Zupp =MiN(Zs, Ve,max /Vops — 1). While zio, is determined by 
the minimal value of ve, i.e., zioy 2max(0, Ve min/Vobs — 1). 
The optimal SNR of the background emission for an obser- 


vation time T is given as 
9HT qe 

Sm =| S 

0 


1/2 


Yrs) (vob) 
dV 0s 


vos Al (Yobs )S n2(Vobs) 
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where S j1(Vobs), S 12(voos) are the noise power spectral densi- 
ties of the two detectors, and y is the normalized overlap re- 
duction function. For two colocated and coaligned detectors, 
y = 1, and we simply assume S p1(Vobs) = S jo(Vops). 

By considering different MMFR models (MFRI and 
MFR2) and EOSs (BSk21 and CDDM2), in Fig. we plot 
the SGWBs contributed by an ensemble of newly born mas- 
sive magnetars with Mmg = 2.4743Mg. As mentioned be- 
fore, depending on the EOS, the radii, representative elliptici- 
ties, and initial spin frequencies of magnetars are taken to be 
R = 12.66 (16.31) km, eg = 0.002 (0.005), and v; = 1660.85 
(888.97) Hz for EOS BSk21 (CDDM2), while the dipole mag- 
netic fields are taken the same as By = 10! G for all magne- 
tars [34]. Consequently, the maximal observed GW frequency 
is vps = 3322 Hz for EOS BSk21 and 1778 Hz for EOS 
CDDM2. For the same EOS, the background spectra calcu- 
lated by using MFR1 have different shapes in comparison with 
those derived by using MFR2. For instance, the peak frequen- 
cies are at 830—1690 Hz for MFRI versus 1132—2455 Hz for 
MFR2. Moreover, assuming EOS CDDM2, the spectrum cal- 
culated based on MFRI with pscrg(0) = 1.1 x 1076 events 


3 —L- ~ — — 
10° £—— MFR2 with p:(0)=10 Mpe Myr", CDDM2 E PIE 
[— — MFR2 with 62(0)=10 Mpc“ Myr", BSk21 E p: 4 
af MFR2 with 62(0)=1 Mpe*Myr", CDDM2 ^ poa 
10 Fos. MER2 with 5 (0)-1 Mpc^Myr", BSk21. ai io 1 
M cie M \ 
10°F one 1 1 
E Mg | i 
x y 
$ c" M p j 
(el ede cic at = y ‘4 
10°F l 1 
i 
F od ad 3 
—— MFR1 with p,,,.,(0)=1.1E-6 Mpc" yr’, CDDM2 
10" E — — MFR1 with p,...(0)-1.1E-6 Mpc^yr", BSk21 — 3 
E —— MFR1 with p,,,,(0)=1.15E-7 Mpc°yr', CDDM2 1 
" ---- MFR1 with p... (0)71.15E-7 Mpc yr", BSk21 
10 " L 1 1 1 
10 100 1000 


Vobs (HZ) 


FIG. 3: Dimensionless GW energy density Qgw versus the obser- 
vational frequency vops, calculated for different MMFR models with 
various local formation rates and different EOSs as shown in the leg- 
end. The newly born massive magnetars have Mmg = 2.4743Mo, and 
By = 10P G. Other quantities for the magnetars are taken to be EOS 
dependent (see the text). For comparison, the detection thresholds of 
aLIGO (black dotted line) and the future ET (black dash-dotted line) 
calculated following Eq. (136) in and assuming a 1 yr observa- 
tion time are also shown. 


Mpc ^? yr! (red thick solid line) dominates the spectrum ob- 
tained by using MFR2 with ò¢(0) = 10 Mpc ?Myr ! (green 
thick solid line) at vobs x 1200 Hz, however, succumbs to the 
later above 1200 Hz. This is because MFRI predicts higher 
pwrR(2) at high z, but lower pmrr(z) at low z. The source for- 
mation rate at high (low) z mainly contributes to background 
emission at low (high) frequency band [58]. 

It is obvious that the background spectra are strongly de- 
pendent on the assumed EOSs as shown in Fig. [3] The EOS 
BSk21 leads to a cutoff at ve ~ 2453 Hz in the GW signal 
emitted by a single magnetar with Mmg = 2.4743 Mg (see Fig. 
[D. Such cutoffs also appear in the background spectra emitted 
by an ensemble of such magnetars if EOS BSk21 is assumed. 
However, the observed cutoff frequency depends on the maxi- 
mal redshift of the MMFR model as veut = 2453/(z.-1) = 350 
Hz with z, = 6 for the two MMFR models. In contrast, us- 
ing a commonly assumed mass 1.4Mọ for newly born mag- 
netars, the resultant background spectrum extends down to 
several hertz without a cutoff because these magnetars do not 
collapse [30/33]. With EOS BSk21, the background emis- 
sions cover the frequency band from 2453 to 350 Hz, which 
are not in the sensitive band of ET. The background spec- 
trum may only be detected by ET in the case of MFRI with 
PsanB(0) = 1.1 x 1076 events Mpc^?yr^! (red thick dashed 
line in Fig. B). The corresponding SNR of the background 
spectrum is (S/N)g = 3.38 (red hollow star in Fig. |4) for ET 
when an observation time T — 1 yr is assumed. The SNR is 
slightly above the detection threshold (S/N)g tn = 2.56 of ET 
(dotted line in Fig. Bih. A lower SNR with (S/N)g = 0.69 
(green hollow star in Fig. is obtained using MFR2 with 
pe(0) = 10 Mpc ?Myr !. Consequently, it should be hard to 
confirm EOS BSk21 through direct observation of such a cut- 
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FIG. 4: The SNR (S/N)s of SGWB from a massive magnetar popu- 
lation versus the magnetars’ interior toroidal magnetic fields B,. The 
SNRs are calculated by assuming different EOSs and various MM- 
FRs as depicted in the legend. All SNRs are derived with respect to 
ET for an observation time of 1 yr. The dotted line shows the de- 
tection threshold of ET. The stars represent the SNRs calculated by 
adopting various representative ellipticities (corresponding to various 
EOSs) and for different MMFRs (see the text). 


off in the background spectrum emitted by massive magnetars. 


As a comparison, if the EOS CDDM2 is assumed, the 
resultant background spectra emitted by all magnetars with 
Mmg = 2.4743Mg show cutoffs at voy ~ 124 Hz. The spectra 
extend from 1778 Hz down to 124 Hz. Using ET and taking 
T = 1 yr, the SNR of the background spectrum is 50.49 (red 
filled star in Fig. for MFR1 with pscrg(0) = 1.1 x 1076 
events Mpc?yr !. The SNR is reduced by about an order 
of magnitude for the minimal local SGRB formation rate 
PsGRBmin(O) = 1.15 x 1077 events Mpc^?yr^!. Moreover, 
adopting MFR2 with o?(0) = 10 Mpc ?Myr“!, the resultant 
SNR is 10.28 for ET (green filled star in Fig. [4). The SNRs 
are all above the detection threshold of ET, which suggest that 
the spectra may be detectable by the proposed ET. If the cut- 
off at ~ 100 Hz in the SGWB from massive magnetars could 
be detected in the future, then the EOS of dense matter may 
be consistent with QS EOS CDDM2. Furthermore, the EOSs 
(e.g., BSk21 and APR) which provide relatively small Mroy 
could be excluded. Of course, the nondetection of the back- 
ground emission may have the following reasons: (i) the ac- 
tual pwrn(z) is much smaller than those we adopted here; (ii) 
the EOSs (e.g., BSk21 and APR) which provide relatively 
small Mroy are favorable; (iii) the magnetars actually have 
much lower B+. 

Related to the last point of the above reasons, the effect of 
B, on the detection of the SGWB is shown in Fig. [4] All the 
SNRs of the background emissions are calculated by adopting 
the maximal local event rate of each MMFR model, an obser- 
vation time T = 1 yr, and with respect to ET, which has the 
best designed sensitivity. In all cases, as B, increases, the SNR 
of the background first rises rapidly then slowly. Specially, the 
same as in Fig. |2] when EOS BSK21 is assumed the SNR of 


the background tends to be saturated when B, becomes large 
enough. In this case, even MFRI with psagp(0) = 1.1 x 1076 
events Mpc ?yr-! and the maximal toroidal fields B, ~ 10" 
G are taken, (S/N)g is slightly above detection threshold of 
ET. This just reflects that the background emission from mas- 
sive magnetars is difficult to be detected if EOSs with small 
Mmqov are preferred as discussed above. In contrast, for the 
same MMER and Bj, (S/N)g derived based on EOS CDDM2 
is at least ~ 10 times higher than that obtained by assuming 
EOS BSk21 (see Fig. [4). Hence, if EOS CDDM2 rather than 
EOS BSK21 is preferred, detection of the background emis- 
sion from massive magnetars may be promising. 

The SGWB produced by an ensemble of massive magne- 
tars is expected to be continuous if the duty cycle DC = 


T T(l + Z)omrr(z)dV > 1 is satisfied [59]. The quantity T 
is the duration of the GW signal emitted by a single magnetar. 
The expression for the comoving volume element dV can be 
found in (58]. In the case of EOS BSK21, since the 2.4743M, 
massive magnetar has a dipole field B4 = 10? G and repre- 
sentative ellipticity eg — 0.002, its lifetime can be determined 
to be r = 140.50 s. The massive magnetar has a much longer 
lifetime of r = 995.86 s if EOS CDDM2 is assumed, even 
though its representative ellipticity (eg = 0.005) is larger in 
this case. Using MFR2 with p?(0) = 1 Mpc?Myr ! (the low- 
est MMFR), we have DC ~ 39.02 (276.56) for EOS BSk21 
(CDDM2), which means that the SGWB from these mag- 
netars is continuous. We note that even when EOS BSk21 
and ultrastrong toroidal fields B; ~ 10!7 G are assumed for 
the massive magnetars, the produced SGWB is still contin- 
uous if the MMFR is not as low as MFR2 with 62(0) = 1 
Mpc ^? Myr !. 


V. CONCLUSION AND DISCUSSIONS 


As one of the most promising targets for GW detection, 
newly born magnetars, if produced by NS-NS mergers, their 
masses should be much larger than the generally assumed 
value 1.4Mọ for NSs. The masses of the newly born mag- 
netars, the EOSs, and the MMFR models all have an im- 
pact on the SGWB produced by the massive magnetar pop- 
ulation. By taking into account these effects, we estimated the 
SGWB produced by the newly born massive magnetars. For 
the NS EOS BSk21 and QS EOS CDDM2 adopted here, the 
resultant background spectra contributed by massive magne- 
tars with Mmg = 2.4743Mg show cutoffs at about 350 Hz, 
and 124 Hz, respectively. The frequency ranges of back- 
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ground emissions are different for the two EOSs. Assuming 
EOS CDDM2 and representative ellipticities eg = 0.005 for 
the 2.4743Mg massive magnetars, even using MFRI with the 
minimal local rate, the SGWB contributed by an ensemble 
of such magnetars may be detected by the future ET. While 
using MFR2, the background emission from these magnetars 
may be detected by ET only if the local merger rate satisfies 
pe(0) x 10 Mpc ?Myr'!. However, if EOS BSk21 (and rep- 
resentative ellipticities eg — 0.002) is assumed, the SGWB 
may be detected by ET only when MFRI with the maximal 
local rate is adopted. For the same MMFR and B,, adopting 
EOS CDDM2, the SNR of the background is at least ~ 10 
times higher than that obtained based on EOS BSk21. This, 
in turn, may indicate that if such background emission could 
be detected, EOS CDDM2 should be more favorable. The rel- 
atively low SNR of the background emission indicates that it 
may be unlikely to test EOS BSk21 via detecting the cutoff in 
the background spectrum. However, if the cutoff at 100 Hz 
in the SGWB from massive magnetars could be detected in the 
future, the QS EOS CDDM2 seems to be favorable. In addi- 
tion, successful detection of background emission at = 100 Hz 
could reasonably exclude the EOSs which present relatively 
small Mroy (e.g., EOSs BSk21 and APR). Finally, detecting 
the GW emission during the entire formation process (from 
the final binary inspiral process to the magnetar phase) of a 
newly born massive magnetar may still be an effective way to 
probe the EOS of dense matter because one need not estimate 
the MMER, which is actually rather uncertain. 


Improvements are still needed in order to obtain a more 
realistic SGWB produced by an ensemble of massive mag- 
netars. In this paper, we only assume typical mass Mmg = 
2.4743 Mo for the massive magnetars. Actually, the masses of 
these massive magnetars should distribute in a certain range as 
inferred from the mass distribution of Galactic NS-NS bina- 
ries BHA. In future work, we will combine the mass distribu- 
tion of massive magnetars with various NS/QS EOSs to study 
the SGWB produced by the massive magnetars in detail. 
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